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Abstract 



We obtain a class of soliton solutions of the integrable hierarchy which 
has been put forward in a series of works by Z. Qiao. The soliton 
solutions are in the class of real functions approaching constant value 
fast enough at infinity, the so-called 'dark solitons'. 
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>■ 1 Introduction 

ov 

^ . The interest inspired by the Camassa-Holm (CH) equation and its singular 

peakon solutions |1] prompted search for other integrable equations with 
similar properties. An integrable peakon equation with cubic nonlinearities 
has been discovered first by Qiao [UJ and studied further e.g. in [TJl ITB] . 
Another equation with cubic nonlinearities has been found by V. Novikov [5]. 
The Lax pair for the Novikov's equation is given in [7], (see also a remark on 
the peakons of Qiao's equation in [?]). Actually the Qiao's equation 

m t + (m(u 2 - u 2 x )) x = 0, m = u-u xx (1) 

together with the CH equation 

m t + 2u x m + um x = 0, m = u — u xx (2) 

belong to the bi-Hamiltonian hierarchy of equations described by Fokas and 
Fuchssteiner [3]. The Qiao's equation has a distinctive W/M-sh&pe travelling 
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wave solutions [HI [12]- The peakons of Novikov's equation have been studied 
in [8]. 2 + 1 dimensional generalizations of Qiao's hierarchy are studied in 
[3]. Single peakon, mutil-peakon dynamics, weak kink, kink-peakon, and 
stability analysis of the Qiao's equation were studied in [H] and [5]. For the 
CH and related equations one can consult the monographs [6j EJ [10] and the 
references therein. 

Equation ([T]) can also be written as 

m t + (u — u x )m x + 2u x m = 0. (3) 

Qiao presented a 2 x 2 Lax pair for this equation given by the linear 
system & x = U*, * 4 = V* with 



U 
V 



— \ \m\ 
\m\ \ 1 ' 



A~ 2 + \{u 2 - ul) -\-\u - u x ) - \m\(u 2 - u 2 x 

\^{u + u x ) + \m\(u 2 - ul) -\- 2 - \{u 2 - u 2 x ) 

There is another equation from the same hierarchy, 



• ( —) -(—,) =0- (■>) 



m 2 J „ V m 2 



for which the V-operator is Sf^ = where 



(4) 



\2 I m(in x -m xx )+3ml 
V„ — — I m A m 4 1 (f.) 

2 o I ,7 m(m x +m xx )-3ml A | • V u 7 
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The (white) soliton solutions of ([T]) and §5§ have been found previously 
[T5| IT7] . These studies rely on the fact that the spectral problem for ([1]) is 
gauge-equivalent to the one for the mKdV equation. In this study we will 
present soliton solutions approaching a constant value for \x\ — > oo (dark 
solitons). To this end we are going to formulate the spectral problem in the 
form of a Schrodinger operator, which is the same spectral problem as for 
the KdV equation. 

2 Reformulation of the spectral problem 

Let us consider solutions such as 

m(x,t) > 0, lim m(x,t) = m , (7) 

21— »±00 



2 



where tuq is a positive constant. Let us assume also that m(x, ■) — m e S( 
for any value of t. One can reformulate the spectral problem into a scalar 
one as follows. Introducing ^ = (■0, 0) T the matrix Lax pair written in 
components is 



2^ 



With a change of coordinates 



9 = —d 

w y w x ) 

m 



—mXip + 0. 



(8) 



we obtain the following scalar spectral problem for (f>(y, A) (sometimes we do 
not write the argument t which is an external parameter for the considered 
spectral problem) 



+ 



= A 2 0- 



(9) 



Note that this is a Schrodinger's operator with a potential 



U(y,t) 



(10) 



It is well known how to recover U(y,t) from the scattering data of (Q, how- 
ever the solution is m(y,t) and its recovery from U(y,t) necessitates solving 
a nonlinear (Riccati) equation. We can express m(y,t) in terms of the eigen- 
functions of the Schrodinger's operator. We introduce p(y, A) = ^ from 
which we immediately obtain 



Py + P 



U(y) - A 2 



If we define 
then we have 



Po(y) = p(y,o) 



U{V) = P0,y + Po- 

However, due to (fTUl) we now have — = p or 



m(y,t) 



1 



<f>{y,t,\) 



Pa(y,t) (j) y (y,t,\) 



A=0 
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So far we treated y as a new variable instead of x. However we can treat 
y as a parameter, and then ffTTT) represents the solution in parametric form, 
where the original variable x is given due to ([8]), (TTTj) by: 



x(y,t) = 2\n4>(y,t,0) + const. (12) 

Assuming that <p(y, t, 0) is everywhere positive, we have a solution in para- 
metric form ( II ip . (j!2p given entirely in terms of the eigenfunctions 0(?/, i, 0). 
One can write formally the solution (neglecting the constant in the last for- 
mula) as 



m(x, t) 



6(x-2ln<P(y,t,0))dy. 



(13) 



3 Inverse scattering and Soliton solutions 

From (j7j) and (ITU1) it follows that U(y) does not decay to when y — > ±oo. 
To this end we introduce the modified potential 

1 



U(y) = U(y) 
for which lini| 2/ |_ H30 U(y) = 0. So we have 



m 



2 • 



'iW 



U{y) - 



mi 



or, introducing a new spectral parameter 



A 2 - A 



we have a standard spectral problem 



- <p yy + U(y)<p(k, y) = k 2 4>(k, y), U{y) e 5(R). 



(14) 



(15) 



(16) 



When A = however we find = ±— for k. This means that if one takes 

mo 

an eigenfunction 0(fc, y) of ( fl6|) analytic in the upper (lower) half complex 
/c-plane, one should evaluate it at k — — (k = — — ): 

1 ' ran \ mn ' 



m(y,t) 
x(y,t) 



(f>(y,t,k) 



<f>v(.V,t, k) 



(17) 



k=±- 



21n(/ 
4 



y,t,± — 

m 



The spectral theory for the problem ffTBj) is well developed, e.g. [IB]. We 
are going to use these results to construct the soliton solutions of ([T]), fl5j). 
One can introduce scattering data as usual. For the time-dependence of the 
scattering data one needs the time-evolution of the eigenfunction (f>(k,x). 
The Lax-pair in x and t variables for ([I]) has the form 



m x 




m x 




m 




2m 


1 


u x -\- u xx 




u + u x 




m 


<f>- 


X 2 m 



+ 



u 



u 



21 



<Px + 70, 



(19) 



(20) 



where 7 is an arbitrary constant. The second equation, (120]) in asymptotic 
form x — > ±00 is 



->■ - 



1 

A 2 " 



<Px + 70, 



or, in terms of k, y-variables when y — > ±00, 



fc 2 mg + 3 
/c 2 mo + 1_ 



0y + 70, 



(21) 



since 



lim m 

|j/|-k» 



lim u 

|l/|->oc 



mo- 



Defining Jost solutions by 



lim ^(y^e^l, 

y— >±oo 

such that 

V-(y> fe ) = a (. k )v+(y, k ) + b(k)<p+(y, k), fceK 
and noting that </?_ — >■ ae~ lky + 6e ifc?/ when y — >■ 00 we find from (1211) 



rrin 



a t 



b t = - 

Requiring a t = 0, we find 

b t = 



k 2 m 2 Q + 3 
k 2 ml + 1 
k 2 ml + 3 
/c 2 7Bq + 1 

k 2 m\ + 3 
k 2 m 2 n + 1 



zfca) + 7a 
zfcfr) + 76. 



(22) 
(23) 

(24) 
(25) 



-2K — - 

2 



6(fc,i) 
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and thus for the scattering coefficient r = b/a we have 
r(k, t) = r(k, 0) exp 



w,q f k 2 ml + 3 



2 \k 2 m 2 + l 

and for the analogue on the discrete spectrum k = in n , 



(26) 



Rn(t) 



b{iK r . 



xa [iK r 



i? n (0)exp 



K n ml(3 - K 2 n ml) 
2(1 - K 2 n ml) 



(27) 



For the equation (J5J) the time evolution of the spectral eigenf unctions is 
given by 



v 2 

<f>t = —4> + (A 2 - v) —<p x + 70, 
m m 

and analogous considerations give 



m(m x + m xx ) - 3m 2 



r(k,t) = r(fc, 0)exp 
R n (t) = i?„(0)exp 



rrr 



-2ik { + k' A ] t 
mi 



1k 7 



m. 



2 K n 1 t 



(28) 

(29) 
(30) 



It is convenient to introduce a dispersion law for the hierarchy, which for 
the considered two members is 



k?t!,q(3— K 2 m.g) 
f( K ) = I 2( 1 -« 2m o) 

^ -m 2 K 2 ) foreq. ©. 



for eq. ([T]) 



Then for the whole hierarchy we can write in general 

R n {t) = R n {0)exp{f{K n )t). 
For further convenience we introduce 



£n = y- 



-t 



2/w 



In 



Rn(0) 



(31) 



The eigenfunctions of the spectral problem (1T6"1) are well known, see e.g. 
[T6] . In the purely iV-soliton case the eigenfunction, analytic in the lower 
complex &;-plane is the Jost solution ip + (y,k) defined in (T22l) which has the 
form 



6 



^V,t,AO = e^(l + X;^) (32) 

n=l 



with the residues T n (y,t) satisfying a linear system 



r n ( y ,t) = iR n {t)e-^y U+ijr, 

\ m=l 

The time-dependence of the scattering data is given by (jHTl) . The iV- soliton 
solution then is given in parametric form by (JTTJ) and (ITS]) for the eigenfunc- 
tion ( |32|) . The condition < K n < Wq 1 is sufficient to ensure smoothness of 
the solitons. 



4 Example: One-Soliton Solution 

The one-soliton solution corresponds to one discrete eigenvalue ki = m\, 
where K\ is real, positive and K\ < m^ 1 . The eigenfunction in this case is 
©2} 

v + (v, t, k) = (i + ■ m iX 2 T ) ■ ^ 

y k — iK\ l + 2 lW e~ 2 " 1!l J 

Evaluated at k = — we find 

mo 

-i, jl A 1 i?i(t)e- 2,tlS/ 
V+Cl/,*, — ) = e m " 1 - — 



From (ITT)) and (ITS)) we obtain the one-soliton solutions 

,(»,«) = ^+2h(l- — ^Sf^ Y (34) 
m \ {1 + Kim ) cosh Ki^i J 

m (y^ = ; ■ (35) 

1— rrtoKi tanh reifi 



The extremum (minimum) of m occurs when 

6 = — lnp-^V 

4ki \l + m Ki/ 



This is a constant value, e.g. the soliton moves with a velocity that 
depends on the dispersion law (i.e. the chosen equation from the hierarchy). 
The profile of the dark soliton is given on Fig. UJ 
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1.85- 



1.80- 



1.75- 




Figure 1: One soliton profile, mo = 2, K\ — 0.2. 

5 Example: Two soliton solution 

In the case of two discrete eigenvalues we compute 



<p+(y,t, — ) = e™o- 

m 



(36) 



where the following notation is utilized: 



y. - HiS 1 

mo J 



J = 1,2. 



From (|17p and (jTSJ) we obtain the two-soliton solutions: 
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Figure 2: Snapshots of the two (dark) soliton solution of the Qiao equation 
([1]), for three values of t : —30, —12 and 30. The other parameters are 
m = 2, Ki — 0.1, k 2 = 0.25. 



x{y,t) = — + 21n^i (37) 
m A 2 

™(y,t) = t^aT- (38) 



where the following notations are used: 

A 1 (y, t) = 1 + e" 2 ^ 1 + e" 2 ^ 2 + ( ^l^l) * e -2«Ci-2«2& 

A 2 (y, t) = 1 + i^T 2 ** + z/ 2 e" 2 ^ 2 + f *^J*) ' ^ ^-2^1-2^ 

m + «i m + k 2 

+ 8(K! - /t 2 ) 2 e -2*ifa-2*afe 

molmo" 1 + K i)( m o 1 + K s) 



7B 1 + K 2 V K 1 + K 2 



+ 4 ^ 2 f*^) e-^~^. (39) 
The interaction of two dark solitons is illustrated on Fig. |2J 
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6 Conclusions 



In this paper we demonstrated how the spectral problem for the Qiao's hier- 
archy can be reduced to the one for the standard Schrodinger operator and 
hence the soliton solutions ('dark' solitons) can be obtained in a straight- 
forward manner. This necessitates constant boundary conditions for the 
solution and also a restriction on the discrete eigenvalues < n n < rn^ 1 . 
It is interesting what happens to the solutions if this condition is violated. 
Based on the similarity with Camassa-Holm equation it is likely that there 
are breaking waves present in this case. Moreover, the equation (Q} has a 
conservation law in the form 

X x (x, t)m(X, t) = m(x, 0) 

where X is the solution of 

X t (x, t) = u 2 (X, t) - u 2 x (X, t), X(x, 0) = x. 

It is likely that this conservation law will play an essential role in the study 
of the wellposedness, existence and breaking of the solutions. 
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